The correspondence between long-range and short-range spin glasses. 
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We compare the critical behavior of the short-range Ising spin glass with a spin glass with long- 
range interactions which fall off as a power a of the distance. We show that there is a value of a 
of the long-range model for which the critical behavior is very similar to that of the short-range 
model in four dimensions. We also study a value of a for which we find the critical behavior to be 
compatible with that of the three dimensional model, though we have much less precision than in 
the four-dimensional case. 



o 



CZ3 



I 

o 
o 



o 
o 



PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.-fq 



I. INTRODUCTION 

In the theory of systems at their critical point it is in- 
structive to consider a range of dimensions d, since above 
an upper critical dimension, du, the critical behavior be- 
comes quite simple and corresponds to that of mean field 
theory. Hence it is desirable to understand critical behav- 
ior up to, and just above, d = du- For the case of spin 
glasses,^ where much of what we know has come from 
numerical simulations, this has been difficult because (i) 
the value of rf„ is quite large [d^ = 6 as opposed to 4 
for conventional systems like ferromagncts) and (ii) slow 
dynamics, coming from the complicated "energy land- 
scape" , prevents equilibration of systems with more than 
of order 10^ spins at and below the transition temper- 
ature Tc. Since the total number of spins V is related 
to the linear size L hy V = L'^, for dimensions around 
du (= 6) it is then not possible to study a range of values 
of L, which, however, is necessary to carry out a finite- 
size scaling^''^ (FSS) analysis. 

It has been proposed^ to try to circumvent this prob- 
lem by using, instead, a one-dimensional spin glass model 
in which the interactions Jij fall off as a power of the dis- 
tance, roughly Jij ~ l/l^^i — Tjl" , since varying a in this 
l-d model seems to be analogous to varying c? in a short- 
range models. In both cases there is a range where there 
is no transition (d less than a lower critical dimension c?;, 
cr greater than a certain value ct;), a range where there is 
a transition with non-mean field exponents {di < d < du, 
ai > a > au for a certain (t„ which turns out to be 2/3), 
and a transition with mean field exponents (d„ < d < oo, 
cr„ > cr > 1/2). The advantage of the l-d model is that 
one can study a large range of linear sizes for the whole 
range of a. Consequently, there have been several subse- 
quent studies^"^^ on these models. 

The question that we tackle here is whether this con- 
nection between long-range models in 1-dimension and 
short range models in a range of dimensions is just a 
vague analogy or whether the connection can be made 
precise in the following sense: for a given d is there a value 
of a such that all the critical exponents of the short-range 



model correspond with those of the long-range model (in 
the sense of Eq. (5) below)? We will denote the value of 
(T in Eq. (5) as a proxy for the dimension d. 

A relation between the long-range (LR) and short- 
range (SR) exponents has been proposed in Ref. 8. We 
reproduce their argument here in a more general formu- 
lation. Consider the singular part of the free energy den- 
sity. For a system in d dimensions it has the scaling form 
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where / is a scaling function, t = {T — Tc)/Tc is the 
reduced temperature, h is the magnetic field (for a spin 
glass it is actually the variance of a random field), u is 
the operator which gives the leading correction to scaling, 
yT is the thermal exponent, y^ is the magnetic exponent, 
and yu (< 0) is the exponent for the leading correction 
to scaling. These exponents can be expressed in terms of 
more commonly used exponents, 
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where v is the correlation length exponent, rj describes 
the power-law decay of correlations at the critical point, 
and w > 0. 

We make a connection between the two models by 
equating the singular part of their free energy densities, 
i.e. 
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In order to compare exponents we need to eliminate the 
different prefactors in front of the scaling functions by 
writing everything in terms of the total number of spins 
V where V = L"^ for SR and V = L ioi LR. Canceling a 
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factor of 1/V on both sides gives 



(4) 



Hence, for each of the exponents, the correspondence be- 
tween the LR and SR vakics is 



yLR(CT) 



ysR{d) 



(5) 



We note that in the mean- field regime, 6<(i<oo,2/3> 
a > 1/2, Eq. (5) holds consistently^ for the thermal, 
magnetic, and correction exponents with 



2cr - 1' 



(mean field regime) , 



(6) 
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0,?/LR = 3 - 2a, VSR = 1/2, z^LR = 
1/(2(7- l),a;sR = (d-6)/2, and wlr = 2-3cr. Further- 
more, the exponents also match to first order in 6 — d 
for the SR model and cr — 2/3 for the LR model. Actu- 
ally, Eq. (5) (at least as applied to the thermal exponent 
v = I/j/t) can be derived for all d and a from a super- 
universality hypothesis. 

In this paper we will investigate whether, for d ~ 3 
and 4, we can find a value of a which satisfies Eq. (5) 
simultaneously for the thermal, magnetic and correction 
to scaling exponents. 

One advantage of long-range systems is that the expo- 
nent 77 is known exactly as was first shown by Fisher et 
al.^^ for ferromagnets. The result for spin glasses is 



2<T - 1. 



so Eq. (5) for the magnetic exponent yn {— {d + 2' 
can be written^ 
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which immediately gives us a value of a which acts as 
a proxy for d provided we know ?7sr('^)- For d — 3, 
the value of tjsr, as well as other exponents, has been 
determined accurately by Hasenbusch et al.^^ and we 
use their values here. In particular, they find 77sr(3) = 
—0.375(10), which, according to Eq. (8), corresponds to 
a proxy value a = 0.896. We shall therefore perform 
simulations for this value of a to see if the other expo- 
nents, j/T and Uu, also match those of the d = 3 results^^ 
according to Eq. (5). 

However, for c? = 4, the values of r/sR and the other 
exponents are not known with great precision, so we carry 
here out a careful study of this model here to determine 
them more accurately. We find 77sr(4) = —0.320(13) 
for which the proxy value of tr, according to Eq. (8), is 
a = 0.790. Wc therefore also study this value of cr so see if 
the other exponents match those of the d ~ A simulations 
according to Eq. (5). 



It is also convenient to note that Eq. (5) for the thermal 
exponent ut (= l/z^) can be written 



i^lr{o-) = di^sR{d) , 



(9) 



and, since ?/„ = —uj, the connection between the correc- 
tion to scaling exponents is 
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To summarize, the main goal of this paper is to see if 
there is a single value of a which simultaneously satisfies 
Eqs. (8), (9), and (10) for d = 3 and (with a different 
value of cr) for d = 4. 

The plan of this paper is as follows. In Sec. II we de- 
scribe the model and the observables we calculate. Sec- 
tion HI discusses the finite-size scaling analysis, while 
Sec. IV describes the details of the simulations. The 
results and analysis arc presented in Sec. V, while our 
conclusions are summarized in Sec. VI. 



II. MODEL AND OBSERVABLES 

Wc consider the Edwards-Anderson spin-glass model 
with Hamiltonian 



H — ^ JijSiSj 

(id) 



(11) 



where the Ising spins Si take values ±1 and the quenched 
interactions are independent random variables, the 
form of which will be different for the different models 
that we study. 

The first model is a nearest neighbor spin glass in four 
dimensions in which the Jij take values ±1 with equal 
probability if i and j arc nearest neighbors, and are 
otherwise, i.e. the probability distribution is 




\ [5{J^j - 1) + 5{Jij + 1)] , neighbors), 

(otherwise). 

(12) 

The advantage of the ±1 interactions is that we are able 
to use multispin coding, in which the interactions and 
the spins are represented by a single bit rather than a 
whole word. In fact, our C code uses 128-bit words, us- 
ing the streaming SIMD extensions, so we simulate 128 
samples in parallel. In order to gain the full speedup, 
we use the same random numbers for each of the 128 
samples in a "batch". Hence, while the results for each 
sample are unbiased, there may be correlations between 
samples in the same batch. Consequently, when we es- 
timate error bars we first average over the samples in a 
batch and use this average as a single data point in the 
analysis. Data from different batches are uncorrelated. 

The spins are on a 4-dimensional hypercubic lattice 
of linear size L with periodic boundary conditions. The 
total number of spins is V ~ L'^ . 
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The description of the interactions wc take for the 1-d 
models is a bit more comphcated. The interactions must 
fall off with distance such that 

oc , (13) 

where r-y = \ri — rj\ (the i = 0,l,...L — 1 sites in the 
graph are placed in a circle of radius L/{2t:), site i is at 
angle i27r/L). On the other hand [■ • • ]av denotes an aver- 
age over the interactions. The simplest way to do this is 
to have every spin interact with every other spin with an 
interaction strength which has zero mean and standard 
deviation oc l/rfj. However, this is inefficient to simu- 
late for large sizes, because the CPU time per sweep is 
of order L^, rather than Lz in short-range systems with 
coordination number z. Fortunately, it was realized by 
Lcuzzi et al.,^° that one can have the CPU time scale 
also like Lz for the long range model if one dilutes it. In 
their version, most interactions are zero and those that 
are non-zero have a strength of unity (i.e. the strength 
does not decrease with distance). Rather it is the prob- 
ability of the interaction being non-zero which deceases 
with distance. In the specific construction of Leuzzi et 
al.^^ there are a total of Lz/2 non-zero interactions with 
an average degree (i.e. coordination number) of z and the 
probability of a non-zero interaction given by 

= 1 - exp(-A/r2/) (~ A/rl^ at large r,,) , (14) 

where A is chosen so that the mean degree is equal to 
some specified value z. 

In the Leuzzi et al model, the degree is not the same 
for all sites but has a Poisson distribution with mean 
z. Since we wish to implement multispin coding, and 
since the computer code for this depends strongly on the 
degree (and gets complicated for large degree) , we study, 
instead, a model with fixed degree. 

We are not aware of any simple algorithm to generate 
bonds of arbitrary length such that each site has a spec- 
ified number of bonds (z here) and the probability of a 
bond between i and j varies with distance in some 
specified way (oc 1/r^j' here). We therefore construct the 
Hamiltonian for which we will simulate the spins by first 
performing a Monte Carlo simulation of the bonds. A 
similar (but simpler) problem was resolved in this way in 
Ref. 20. We take the "Hamiltonian" of the bonds to be 
given by 

k I 

where eij = or 1, in which 1 represents a bond present 
between sites i and j, and represents no bond. Graph- 
ically, we regard each site i as having z "legs" associated 
with it, and we initially pair up the legs in a random 
way, representing each connected pair graphically as an 
"edge" and giving the value etj = 1 to all edges while 
all other pairs {i,j) have Eij = 0. We then run a Monte 
Carlo simulation in which the non-zero are swapped 
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FIG. 1: Each site has a fixed number of "legs" (here we show 
three) and these legs are paired up by "edges". In the top 
row, one edge connects sites 1 and 2, and another edge con- 
nects sites 3 and 4. A basic Monte Carlo move for the bond- 
generation simulation consists of reconnecting two edges, as 
shown in the bottom row. (Other edges are present but not 
shown.) 



according to a Metropolis probability for the Hamilto- 
nian in Eq. (15). To maintain exactly z non-zero e's for 
each site the basic move involves reconnecting two bonds 
as shown in the sketch in Fig. 1. 

Specifically, we first choose site 1 in Fig. 1, with uni- 
form probability among the L possible choices. Next, site 
3 is chosen with probability proportional to l/r^"" (r is 
the distance among sites 1 and 3). Finally, site 2 (site 4) 
is chosen with uniform probability among the z "neigh- 
bors" of site 1 (site 3). Before the move is attempted, 
wc need to check that the sites 1, 2, 3 and 4 verify two 
consistency conditions. First, the four sites should be all 
different. Second, we require that neither sites 1 and 4, 
nor 2 and 3, are paired. If the consistency conditions are 
met, the basic move can be attempted and then be ac- 
cepted/rejected with Metropolis probability. One sweep 
corresponds to Lz selection of sites of type "1" in Fig. 1. 

After a suitable equilibration time,'^^ we freeze the e^j, 
and the resulting set of non-zero e^j defines a "graph" . 
Each of the 128 samples in a single batch of the multispin 
coding algorithm has the same graph. On the edges of 
the graph wc put interactions with values ±1 with equal 
probability chosen independently for each edge in each 
sample in a batch. The result is that the probability 
distribution for a single bond is given by 

P(J,,) = (l-p.,)<5(J,,)+p.,i [5{J,j-l) + SiI,, + l)] , 

(16) 

in which pij is given approximately by Eq. (14) for an 
appropriate choice of A corresponding to the specified 
value of z. However, the bonds are no longer statistically 
independent; rather there are correlations which ensure 
that each site has exactly z non-zero bonds. For both 
a = 0.896 and 0.790 we take z = 6 neighbors. 

We now describe the quantities that we calculate in 
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the simulations. The spin glass order parameter is 



(17) 



where "(1)" and "(2)" are two identical copies of the 
system with the same interactions. Its Fourier transform 
to wavevector k is denote by q{k). We will calculate the 
spin glass susceptibility 

XSG = ^[(<7')]av, (18) 

and also its wavevector-dependent generalization, 

XSG(fc) = ^[(k(fc)P>]av. (19) 

From this we can extract the correlation length, ■^'^^"^■^ 
1 



/ X(0) 

2sin(7r/L) V x(fci) 
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(20) 



where fei is the smallest non-zero wavevector, ki = 
(27r/L)(l, 0,0,0) for the 4-d model and fci = 27r/i for 
the long-range models in 1-d. Other quantities that we 
calculate, important because they are dimensionless like 
£,l/L^ are the moment ratios. 



C/4 



2\12 
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2\2 ■ 



and the susceptibility ratio 
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av 



Ri 



Xsajki) 

XSG{k2) ' 



(21) 



(22) 



(23) 



where k2 is the second smallest non-zero wavevector, 
k2 = (27r/L)(l, 1,0,0) for the A-d model and ^2 = A-k / L 
for the long-range models. We will also determine deriva- 
tives with respect to /? of several of these quantities using 
the result 



(O) {-H). 



(24) 



III. FINITE-SIZE SCALING ANALYSIS 

Using data from finite-sizes, we have to extract the 
transition temperature Tc, the correction to scaling ex- 
ponent w (since corrections to scaling are significant), the 
correlation length exponent and (for the short-range 
model which it's value is not known exactly) the expo- 
nent 77. In this section we show how to include the leading 
correction to FSS. There are several sources of subleading 
corrections which will not be included in the formulae in 



this section, though we will try to include them empiri- 
cally in some of the fits to the data, as discussed later in 
the section. 

It is desirable to compute the various quantities one 
at a time so the value of the exponents depend on each 
other to the least extent possible. We therefore adopt 
the following procedure. 

We start with the finite-size scaling (FSS) form of a 
dimensionless quantity, since these quantities are simpler 
to analyze than those with dimensions and so they form 
the core of our analysis. 

Dimensionless quantities are scale-invariant, which 
means that at Tc they remain finite (neither zero nor infi- 
nite) in the limit of large L. However dimensionless quan- 
tities are not only scale-invariant, they are also universal 
(i.e. they remain constant under Renormalization- Group 
transformations). Examples of dimensionless quantities 
are 1/4^ U22 and i?i2. The distinction among scale- 

invariant and dimensionless quantities has been stressed 
in Ref. 18. Here we will discuss dimensionless quantities, 
but will comment on quantities which are scale-invariant 
but not dimensionless in the last paragraph of this sec- 
tion. 

A dimensionless quantity f{L,t) has the FSS scaling 

fQj.j^23-25 

f{L,t) ^ Fo{L'/''t) + L-'^F,{L'/-t), (25) 
where tu is the correction to scaling exponent, and 

T-Tr. 



t = 



(26) 



We are interested in the behavior at large L and small t, 
and including just the leading corrections in 1/L and t 
gives 

/(L,t) ~ Fo(0) + Li/'^ti?^(0) + L-"Fi(0). (27) 

It will be useful to determine the values of t*j^ where the 
quantity / takes the same value for sizes L and sL, where 
s is a scale factor which we shall take to be 2 here. We 
have 

Fo(0) + (sLy/'tl F^(0) + (,5L)-"i?i(0), (28) 
which gives 



or equivalently, to leading order. 



■^c - ML 



/3c 



where the non-universal amplitude is given by 

,/ _ (l~.s-")Fi(0) 
^ (sV--l)F^(O)- 



(29) 



(30) 



(31) 
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One can use Eq. (30) to locate /3c. As we shall see, the 
exponents lo and 1/v are determined separately, and we 
use those values when fitting the data to Eq. (30). 

We shall determine the critical exponents using the 
quotient method, which is a more modern form of 
Nightingale's phenomenological renormalization.^^ First 
we determine the correction exponent uj by applying the 
quotient method to dimensionlcss quantities. Consider a 
second dimensionlcss quantity g{L,t) which varies near 
Tc in the same way as / in Eq. (27), i.e. 

g{L,t) ~ Go(0) + L^^-'tG'a(Q) + L-"Gi(0). (32) 

Now compute g{L,t) at t'^, given by Eq. (30), the tem- 
perature where results for L and sL intersect for some 
different dimensionlcss quantity /. We have 



g{L,tl)^Go{0)+Ai^fL- 



(33) 



where A^'f = A(G'o{0) + Gi(0). While this could be 
used directly to determine lo it is more convenient to take 
the ratio (quotient) of this result with the corresponding 
result for size sL, i.e. 



= 1 + BI-fL- 



(34) 



where the amplitude Sf'^ is non-universal (because of 
the definition, it is zero if the quantities / and g are the 
same). Eq. (34) is the most convenient expression from 
which to determine u since it just involves the one un- 
known exponent w, and one amplitude B. These quanti- 
ties can be determined by a straight-line fit to a log-log 
plot oi Q{g) — 1 against L. 

To determine the other exponents v and 77 we need to 
consider the FSS scaling form of quantities which have 
dimensions. Consider some quantity O which diverges in 
the bulk like t~^° . Including the leading correction it 
has the FSS form 



-l/u. 



(35) 



where yo = xq/v. Repeating the above arguments, and 
determining O for sizes L and sL at the intersection tem- 
perature t*j^ for the dimensionlcss quantity / for sizes L 
and sL, the quotient can be written as 



B?'fL- 



(36) 



Using the value of a; determined from Eq. (34) the expo- 
nent yo is determined from Eq. (36) by a straight line fit 
to a plot of Q{0) against 1/L'^. 

To determine 77 we can use Eq. (36) for the spin-glass 
susceptibility xsGi since yo = — rj because the suscep- 
tibility exponent 7 (= x^^^) = (2 — ri)u. To determine u 
we note that £,l/L is dimensionlcss and so has the same 
FSS scaling form as in Eq. (25). Differentiating, for in- 
stance, S,L with respect to /3 brings down a factor of L^^'^ 



and so yo = 1 -I- l/i^ in this case (yo = \/v li we take 
the logarithmic derivative). Hence we determine l + l/v 
from Eq. (36) with O given by the /3 derivative of . 

To conclude, to carry out the FSS analysis we do the 
following steps: 

1. Determine w from Eq. (34) for one or more dimen- 
sionlcss quantities /. 

2. Using the value of uj so determined, obtain 1 + 
(and 2 — 7] where necessary) from Eq. (36) with 
O = xsG a-nd O = 9^^/9/3 respectively. 

3. Using the value of w from stage 1 and l/iy from 
stage 2, determine /3c from Eq. (30). 

The error bars for l + l/i' and 2 — 77 from stage 2 will have 
a systematic component, coming from the uncertainty 
in the value of w from stage 1, as well as a component 
from statistical errors in the data being fitted. Similarly 
the error bar in /3c from stage 3 will have a systematic 
component due to uncertainty in the value of a; -I- 1/v. 

Each of these three stages only requires a straight-line 
fit. However, in practice things are a little more tricky. 
We would like to use data for as many sizes as possi- 
ble, but in practice the smaller sizes arc affected by sub- 
leading corrections to scaling so we can only use data 
for the larger sizes. It is therefore necessary to include 
only a range of sizes for which the quality of the fit is 
satisfactory. 

In some cases we try to incorporate a sub-leading cor- 
rection to scaling to increase the range of sizes that can 
be used. These arc of different types, one of which is 
higher powers of the leading correction, and this is the 
only one we will include here in order to avoid introducing 
too many additional parameters. In other words, when 
we include sub-leading corrections we will do a parabolic, 
rather than linear, fit to the data as a function of 

In order to increase the number of data points rela- 
tive to the number of fit parameters, we will often do a 
combined fit to several data sets. For example, when es- 
timating UJ we will determine the /3£ from one dimension- 
lcss quantity /. and then determine two (or more) other 
dimensionlcss quantities at these temperatures. These 
data sets will be simultaneously fitted to Eq. (34) with 
the same value for a; (since this is universal) but differ- 
ent amplitudes B (since these are non-universal). Hence, 
by combining two data sets, we double the amount of 
data without doubling the number of fit parameters. It 
should be mentioned that, for a given size, the data 
for the different data sets is correlated, and best esti- 
mates of fitting parameters are obtained by including 
these correlations. In other words, if a data point 
is {xi, yi), and the fitting function is u{x), which depends 
on certain fitting parameters, we determine those param- 
eters by minimizing 



X 



J2[y^~H^^)] {c-'),^ [y^-nix,)], (37) 



6 



where 



(38) 



is the covariance matrix. If there are substantial cor- 
relations in many elements, the covariance matrix can 
become singular, but we have checked that this is not 
the case for the quantities we study. 

We end this section by discussing the FSS of a scale- 
invariant (but dimensionfull) quantity, which turns out to 
be useful in our study of the LR model. Take Eq. (35) and 
imagine that we know exactly the exponent yo- Then, 
0{L,t)/ L^" is scale- invariant, since it remains finite at 
t = even in the limit of large L. This is precisely 
the situation in the LR model, if we take for O the 
SG susceptibility, because, as explained in the introduc- 
tion, the anomalous dimension is a known function of 
a for those models. Nonetheless, Eq. (25) needs to be 
modified when applied to xsg/L'^'^~^ , because the mag- 
netic scaling field u{h,t) is not exactly h, as assumed 
in Eq. (1) (see e.g. Refs. 3,18). Rather, there is a 
non linear dependency on the thermodynamic control 
parameters t and h: Ufi{h,t) = huhit) + 0{h^), where 



Uhit) = 1 + Cit 
Eq. (25) reads 



Cat" 



Hence, the analogue of 



^^^^ = ulit) [OoiL^^^t) + L-"Oi(ii/''t)] . (39) 

We note that the multiplicative renormalization uf^{t) 
cancels out when looking for crossing points, namely 



I2a 



(sL) 



2CT-1 



(40) 



so scales as in Eq. (30) . Unfortunately, the multiplica- 
tive renormalization can no longer be ignored when we 
compute 1/j/ from d^xso/L^'^'^. Indeed, differentiating 
Eq. (39) with respect to /3 and neglecting terms of order 
1/1'^+^^", we find 



9pXSG{L,t) 



L 



ul{t)0',{L^/'^t) 
-ul{t)d[{L'/''t) 



(41) 



rather than Eq. (35). Both Uh and u'^ behave as 
L-independent constants (up to corrections of order 
l/L'^+^Z'') when evaluated at the crossing point given 
in Eq. (30). Hence, the quotient of the /? derivative of 
logxsG is given by 



Qidp^ogxso) = s 



- .1/- 



(42) 



instead of Eq. (36), showing that there are corrections of 
order L^^l^ as well as . For some values of cr, and 
also the 3-(i SR model, one finds Xjv < w so the L~^l^ 
correction dominates. 



TABLE L Parameters of the simulations of the 4-ci model: 
is the number of temperatures with /3max the largest and 
/3min the smallest. The number of Metropolis sweeps is given 



by A^swccp , and the number of samples is 


samp ■ 




1^ -/Vswccp 




/3max 


/3min 


-/Vsamp 


4 2.56 X 10"* 


23 


0.5025 


0.4 


2^" 


5 2.56 X 10^ 


23 


0.5025 


0.4 


220 


6 2.56 X 10^ 


23 


0.5025 


0.4 


220 


8 2.56 X 10^ 


23 


0.5025 


0.4 


220 


10 2.56 X 10^ 


23 


0.5025 


0.4 


220 


12 2.56 X 10^ 


23 


0.5025 


0.4 


220 


16 5.12 X 10* 


23 


0.5025 


0.4 


220 



TABLE II: Parameters of the simulations of the 1-d model 
with (7 = 0.790. See Table I for an explanation of the symbols. 



L 


^ sweep 




/3m ax 


/3m in 


samp 


512 


10" 


16 


0.671 


0.538 


64000 


1024 


10® 


16 


0.671 


0.538 


64000 


2048 


10® 


16 


0.671 


0.538 


64000 


4096 


1.28 X 10® 


16 


0.671 


0.538 


64000 


8192 


1.28 X 10® 


16 


0.671 


0.538 


64000 


16384 


2 X 10® 


16 


0.671 


0.538 


64000 


32768 


2 X 10® 


16 


0.671 


0.538 


64000 



IV. SIMULATION DETAILS 

For each size and temperature we simulate four copies 
of the spins with the same interactions. By simulat- 
ing four copies we can calculate, without bias, quantities 
which involve a product of up to four thermal averages, 
such as the spin glass susceptibility, Eq. (18), the C/4 mo- 
ment ratio, (21), and derivatives of these quantities with 
respect to /? calculated from Eq. (24). 

The simulations use parallel tcmpering^^ (PT) to speed 
up equilibration. For the same set of interactions we 
study Nji values of /3 between /3max and /3min- To ob- 
tain good statistics we simulate a large number, A^sampj 
of samples, where iVgamp is a multiple of 128 because 128 
samples are simulated in parallel by multispin coding. 
For the long-range models there are iVsamp/128 different 
graphs, but each sample for the same graph has different 
interactions. We run for A^swcop single-spin flip (Metropo- 
lis) sweeps performing a parallel tempering sweep every 
10 Metropolis sweeps. The parameters used for the dif- 
ferent models are shown in Tables I-III. 



TABLE III: Parameters of the simulations of the 1-d model 
with (7 = 0.896. See Table I for an explanation of the symbols. 



L 


-^swcep 




/3m ax 


/3m in 


-/Vsamp 


512 


1.28 X 10" 


16 


1.5 


0.6 


12800 


1024 


2.56 X 10® 


13 


1.2 


0.6 


12800 


2048 


1.024 X lO'^ 


14 


1.2 


0.65 


12800 


4096 


8.192 X lO'^ 


16 


1.2 


0.65 


12800 


8192 


8.192 X lO'^ 


16 


1.1 


0.71 


12800 



7 




Niygg 0.4 0.42 0.44 0.46 0.48 0.5 

P 



FIG. 2: (Color online) The difference in the value of the 
^l/L between measurements obtained in the range of sweeps 
Nmcs/'2 to A^MCS and measurements in the range A''mcs/4 to 
for values of A^mcs increasing by factors of 2 up to 
A^swcoD ~ 8.192 X 10^. The data is for the long-range model 
0.896 at 13 



* sweep 

with a 



1.2, the lowest temperature studied. 



To check that the simulations were run for long enough 
to ensure equihbration we adopted the following pro- 
cedure. We divide the measurements into bins whose 
size varies logarithmically, the first averages over the 
last half of the sweeps, i.e. between sweeps iVgweep and 
^swcop/2, the second averages between sweeps 7Vsweep/2 
and 7Vsweep/4, the third between sweeps A^swGop/4 and 
etc. We require that the difference between the 
results in the first two bins is zero within the error bars, 
where we get the error bar for the difference by forming 
the difference between the results for the two bins sepa- 
rately for each sample before averaging over samples. In 
most cases, to be on the safe side, we actually require 
that the differences between the first three bins arc all 
zero within errors. 



FIG. 3: (Color online) A global view of the data for the cor- 
relation length divided by L for the 4-d model. 
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0.34 
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0.49 0.495 
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FIG. 4: (Color online) An enlarged view of the data in Fig. 3 
showing the region of the intersections. 



V. RESULTS 



This procedure is illustrated in Fig. 2 which shows data 
for the long-range model with V = 4096, a = 0.896 at 
/3 = 1.2, the largest /3 value that we studied. The vertical 
axis is the difference in ^l/L between the bin containing 
measurements in sweeps A^mcs/2 to iVMcs and the bin for 
sweeps in the interval 7Vmcs/4 to iVMCs/2, for different 
values of A^mcs up to iVgwoep = 8.192 x lO'', the value in 
Table HI. Since the two points for the largest number 
of sweeps are zero within errors, it follows that the first 
three bins all agree. 



A. Four-dimensional short range model 

Figures 3 and 4 show results for ^l/L defined in 
Eq. (20) and Fig. 5 shows results for the dimensionless 
ratio of moments U4 defined in Eq. (21). The resulting 
inverse temperatures where data for sizes L and 2L 
intersect, i.e. where their quotient Q is unity, is shown in 
Table IV. Results are given for both ^l/L and U4. 

To compute the correction to scaling exponent w we 
determine the quotient of /L at the C/4 crossing and 
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Q(^l/L) where Q(U4)=1 ^-b- 




0.25 



FIG. 5: (Color online) An enlarged view of the data for U4 
for the 4-d model showing the region of the intersections. 



TABLE IV: Inverse temperatures where data for sizes L 
and 2L intersect, i.e. where the quotient Q is equal to unity, 
for (,l/L and the ratio of moments C/4, for the 4-d short-range 
model. 



L 


PI where Q{£.l/L) = 1 


PI whore Q{U4) = 1 


4 


0.49113 ± 0.00009 


0.49725 ± 0.00011 


5 


0.49598 ± 0.00007 


0.50001 ± 0.00009 


6 


0.49825 ± 0.00006 


0.50118 ±0.00008 


8 


0.50012 ± 0.00005 


0.50180 ± 0.00006 



vice versa. These quotients are shown in Table V and 
plotted in Fig. 6. Fitting the largest two pairs of sizes 
for each quantity to Eq. (34) for s = 2 with the same 
exponent w gives 

wsr(4) = 1.04(10), xVdof = 0.99/1 . (43) 

It should be mentioned that the lines in Fig. 6 are not 
separate fits to each set of data but are combined fits 
including the whole covariance matrix. 

We have tried also fits including subleading correc- 
tions to scaling. For instance, considering, in addition, 
the quotient of R12, defined in Eq. (23), at the cross- 
ings of £,l/L and C/4, and fitting the three largest sizes 
to 1 + BiL^^ + B2L~^'^ gives a satisfactory fit with 
uj = 1.29(26), x^/dof = 2.26/5. However we prefer the re- 
sult oj = 1.04(10) since it has been obtained using larger 
lattices (L > 6). 

Next we compute rj from the quotients of XSG, de- 
fined in Eq. (18), at the crossings oi ^l/L and C/4, which 
are shown in Table VI and Figures 4 and 5. Assum- 
ing u} = 1.04(10), a linear fit to Eq. (36) with s = 2 
and the same value of yo 2 — ?/) for both quantities 
gives, for the largest two pairs of sizes, Q = 2^~^ = 
4.949(45)[li4], x^/dof = 0.42/1, in which the numbers in 



FIG. 6: (Color online) The quotient of the dimensionless 
quantity £,l/L of the 4-d model at the U4 crossing (squares) 
and the quotient of f/4 at the ^l/L crossing (triangles). The 
straight lines represent the best fit to Eq. (34) using the 
largest two sizes, with the correction to scaling exponent lj as 
an adjustable parameter. 



TABLE V: Quotients of U4 at the crossings of ^l/L, and vice 
versa, for the 4-d short-range model. 



L 


Q{U4) where Q(^l/L) = 1 


Q{^l/L) where Q{U4) = 1 


4 


1.01675 ±0.00020 


1.02835 ± 0.00033 


5 


1.01311 ±0.00020 


1.02230 ± 0.00033 


6 


1.01112 ±0.00020 


1.01886 ± 0.00033 


8 


1.00822 ±0.00020 


1.01397 ± 0.00033 



rectangular brackets, [■•■], correspond to the errors due 
to the uncertainty in the value of uj. This fit is shown in 
Fig. 7 by the dashed lines. 

On the other hand, a quadratic fit to Q{xsg) = 
Q + ± B2L~^^ using the largest three pairs gives 

Q = 5.039(10)[li^], xVdof = 0.076/1, which is also an 
acceptable fit, shown by the solid lines in Fig, 7. 

If we assume the larger value for w discussed above, 
namely w = 1.29(26) we find that only a quadratic 
fit is acceptable, and the value for Q is Q = 
4.962(30)[6],xVdof= 0.011/1, which is intermediate be- 
tween the two previous values of Q. We can summarize 
all the numbers with the value 

Q EE 2^-'' = 4.994(45). (44) 

The central value is shown as the solid horizontal line in 
Fig. 7, and the error bars are indicated by the dotted 
horizontal lines. Equation (44) gives 

'7sr(4) = -0.320(13) . (45) 

To compute 1/ we have used the quotients for the /3- 
derivative of ^ at the crossings of ^l/L. The values for 
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FIG. 7: (Color online) The quotients of XSG of the 4-d model 
at the crossings of ^l/-!/ (triangles) and C/4 (squares) as a 
function of L~", where uj has already been determined, see 
Fig. 6, and is given by Eq. (43). The dashed lines are the 
linear fit, with a common intercept on the y axis, to the two 
largest pairs of sizes, and the solid lines are the quadratic 
fit to the three largest pairs of sizes (again with a common 
intercept). The intercept is equal to 2^^''. The horizontal 
lines indicate the final estimate and error bars for Q given in 
Eq. (44). This leads to the final estimate for 77 in Eq. (45). 



TABLE VI: Quotients of xsa at the crossings of ^l/L and 
?74 for the 4-ci short-range model. 



L 


Q(xsg) where Q(eL/L) = 1 


Qkxsa) where Q{IJ^ = 1 


4 


4.6464 ± 0.0022 


5.0077 ± 0.0045 


5 


4.7477 ± 0.0022 


5.0368 ± 0.0046 


6 


4.8074 ± 0.0022 


5.0547 ± 0.0047 


8 


4.8673 ± 0.0022 


5.0522 ± 0.0047 



each pair are given in Table VII. Taking uj = 1.04(10) 
we obtain, fitting the three largest pairs, to Eq. (36) for 
s = 2, 

Q = = 3.828(9) [8], xVdof = 0.68/1, (46) 

which gives v = 1.068(4)[3]. Combining the errors we get 
our final estimate for v as 

i^sr(4) = 1.068(7) . (47) 



TABLE VII: Quotients of the j3 derivative of at the cross- 
ings of £,l/L for the 4-ci short-range model. 



L 


Q{dp£,L) whore Q(a/i) = 1 


4 


3.9581 ±0.0024 


5 


3.9340 ± 0.0026 


6 


3.9133 ± 0.0025 


8 


3.8936 ± 0.0031 



FIG. 8: (Color online) The quotient of dpS^L of the 4-ci model 
at the crossing of £,l/L. as a function of , where ui has 
already been determined, see Fig, 6, and is given by Eq. (43). 
The solid line is a linear fit to Eq. (36) using the three largest 
pairs of sizes. The intercept is equal to 2^^^^" . The final 
value of V is given in Eq. (47). 

The data and the fit are shown in Fig. 8 

Finally we estimate /3c by fitting the crossing points 
for £,l/L and to Eq. (30), using the previously de- 
termined values w = 1.04(10) and u = 1.068(7). The 
data has already been given in Table IV and is plotted in 
Fig. 9. We obtain a good fit considering only the (6,12) 
and (8,16) pairs: 

/3c = 0.50256(14)[15], xVdof = 0.24/1 . (48) 

This fit is shown by the dashed lines in Fig. 9. 

We have tried to (roughly) take into account higher 
order corrections to scaling adding a quadratic term in 
j^-uj-i/u^ We obtain a good fit with the pairs (5,10), 
(6,12) and (8,16): 

/3c = 0.50195(34) [1], xV^of = 0.30/1 , (49) 

and this is shown by the solid lines in Fig. 9. We can 
therefore safely take the value, 

/3c = 0.5023(6) =^ Tc = 1.9908(24) (d = 4) , (50) 

as our final result. 

Wc end this section by comparing our results with 
previous computations by other authors. Marinari and 
Zuliani'^*' studied the A-d spin glass with binary couplings, 
finding Tc = 2.03(3), v = 1.00(10) and 77 = -0.30(5), in 
good agreement with our more accurate estimates. Jorg 
and Katzgraber^^ studied a different version of the A-d 
spin glass which is expected to belong to the same univer- 
sality class. They found v = 1.02(2) and -q ^ -0.275(25), 
which are two standard deviations from our estimate. 
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FIG. 9: (Color online) Values of PI, the crossing points for 
^l/L and f/4, for the 4-d model, and fits as a function of 
l/L'^'^^^'^ , for the 4-d model. We used the values of uj and u 
previously determined, see Eqs. (43) and (47). The dashed 
lines are the linear fit, according to Eq. (30), with a common 
intercept on the y axis, to the two largest pairs of sizes, and 
the solid lines are the quadratic fit to the three largest pairs 
of sizes (again with a common intercept). The intercept is the 
critical coupling /3c. The green data points are the estimates 
for Pc for the two fits, Eqs. (48) and (49). 



Jorg and Katzgrabcr also considered the leading eorree- 
tions to scaling, but found an extremely large exponent, 
Lj « 2.5. They were aware that such a large uj is un- 
likely to be correct, and they attributed their result to 
the small lattice sizes that they could equilibrate. 



One-dimensional long range model with 
o- = 0.790 



From Eq. (8) and the value ?7sr(4) = -0.320(13) for 
the i-d model given in Eq. (45), we see that a = 0.790 is 
a proxy for the A-d short-range model, at least according 
to the comparison of the exponents rj (or equivalently of 
the magnetic exponents yn, sec Eq. (5)). In this section 
we will see if Eq. (5) is also satisfied for the thermal 
exponents j/t (for which Eq. (5) can be expressed in terms 
of f as shown in Eq. (9)), and the correction to scaling 
exponents w (= —yu)- Since ?7lr is known exactly, 2 — 
^lrCc) = 2cr — I, see Eq. (7), we can include xsg/L^'^^^ 
as another scale invariant quantity to be studied. 

We focus on and xsg/L'^'^^^^ data for which are 
shown in Fig. 10, and the corresponding crossing points 
are given in Table VIII. Our first task is to try to de- 
termine the correction to scaling exponent uj. We fit the 
quotients of ^l/^, t^4, and U22 defined in Eq. (22), at the 
crossmg of xsg/L'^"'^, including all the (L, 2L) pairs. A 
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L= ' 512 
L= 1024 
L= 2048 
L= 4096 
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L= 16384 
L = 32768 




FIG. 10: (Color online) Correlation length in units of the sys- 
tem size (top) and scale-invariant combination of the SG sus- 
ceptibility and the lattice dimension XSG / L^'^^^ (bottom), 
as a function of the inverse temperature P, for the LR-model 
with a — 0.790. For both quantities, the curves for the differ- 
ent L should cross at temperatures that approach the critical 
point when L grows, see Eq. (30). 



TABLE VIII: Inverse temperatures /3|, , where data for sizes L 
and 2L intersect, i.e. where the quotient Q is equal to unity, 
for xsg/L'^"'^ and £,l/L for the LR model with a = 0.790. 

~ PI where Qixsa/L'"'-') = 1 Pt where Q(^l/L) = 1 



512 
1024 
2048 
4096 
8192 
16384 



0.6538 ± 0.0020 
0.6532 ±0.0018 
0.6516 ±0.0014 
0.6498 ± 0.0012 
0.6500 ± 0.0009 
0.6492 ± 0.0008 



0.6665 ± 0.0066 
0.6598 ± 0.0050 
0.6586 ± 0.0038 
0.6545 ± 0.0031 
0.6541 ± 0.0023 
0.6501 ± 0.0019 



straight line fit, shown in Fig. 11, is acceptable: 

= 0.539(9), xVdof 16.7/14, (51) 

and has a probability of 15%. A quadratic fit to 1 ± 
BiL-"^ + gives a better fit: uj = 0.29(-4 ± 

9),x^/dof = 7/11. This is consistent with the value 
0.26(3) expected from the correspondence in Eq. (10) 
and the value of uj for the 4-d model given in Eq. (43). 
We have also tried fits in which uj is fixed to the value 
0.26. A straight line fit using all the data is very poor, 
X^/doi = 1069/15, whereas a quadratic fit works well, 
X^/dof — 7.5/12, and is shown in Fig. 12. 
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FIG. 11: (Color online) The quotients of dimensionless quan- 
tities ^l/L,U4, and U22 for a = 0.790 at the crossing of 
Xsg/L^'^~^ ■ The straight lines represent the best fit to 
Eq. (34) using all the data, with the correction to scaling 
exponent lj as an adjustable parameter. 

Altogether, we see that our data for the quotients 
of scale invariant quantities do not constrain lj pre- 
cisely. Any value in the range 0.25-0.55 can be consid- 
ered acceptable. Fortunately, this includes the value ex- 
pected from the the correspondence with the 4-d model, 
w = 0.26(3). 

To estimate v we consider the {L, 2L) quotients of the 
logarithmic derivative of xsGi^l, and U4 with respect to 
/?, at the crossings of xsg/L^'^^^ . All these quotients 
should tend to 2^^" for L — s- 00. A straight-line fit ac- 
cording to Eq. (36), allowing w as well as the intercept 
Q to vary, is shown in Fig. 13. The result is 

g = 2^/'' = 1.1703(23), xVdof= 14.24/13, (52) 
wlr(0.790) = 0.277(8) , (53) 

which gives 

t'LR(0. 790) = 4.41(19) . (54) 

This is consistent with the result 4.272(20) expected 
from the correspondence with the 4-d model, see Eq. (9), 
and the the 4-d value of j/ given in Eq. (47), i'sr(4) = 
1.068(7). It is surprising that the fits in Fig. 13 gives 
such a good precision for w, better than using quo- 
tients of scale invariant quantities which we showed in 
Figs. 11 and 12. The result uj = 0.277(8) is consis- 
tent with that expected from the A-d correspondence, 
Lj = 0.26(3). We have also tried a quadratic fit, which 
gives Q = 1.1742(58)[22],xVdof = 9.54/11, and a lin- 
ear fit discarding the L = 512 data which gives Q ~ 
1.1683(15)[62],x7dof = 7.56/8 (both of these fits used 
the value for uj obtained from the correspondence with 



FIG. 12: (Color online) The quotients of dimensionless quan- 
tities ^l/L,U4 and U22 at the crossings of XSG / L^"^^ for 
a — 0.790. The lines represent the best quadratic fit as 
function of l/L", using all the data, where ui is fixed at 
0.26 (= 1.04/4), the value expected from the correspondence 
with the 4-d model, for which the value of uj is given in 
Eq. (43). 



the 4-d model, uj = a;sR(4)/4 = 0.26(3)). These results 
are all consistent with Eq. (54) which we therefore take 
as our final estimate for :^lr (0.790). 

However, the alert reader will recall from Sec. Ill that 
the /3-derivative of xsg/^^"^"^ suffers from two types of 
corrections to scaling, one of order L^'^ and the other 
of order L''^/", see Eqs. (42) and (42). The relation- 
ship between LR and SR exponents in Eqs. (9) and 
(10), combined with our numerical results for the d = 4 
SR- model in Sect. VA, suggests that the two correc- 
tions to scaling are very similar for a ~ 0.790 because 
a;sR,(4) ~ l/t'SR(4). This implies that the two correc- 
tions can be lumped together into a single term to a good 
approximation. Indeed, we have succeeded in analyzing 
our numerical data by considering only the scaling correc- 
tions of order L^". Therefore, although we take Eq. (53) 
as our final estimate for a;LR(0.790); we warn that its 
error is probably underestimated, due to the oversimpli- 
fication in the functional form for the scaling corrections. 

By contrast, we shall see in Sect. V C that for a = 0.896 
the corrections of order L~^l^ turn out to be dominant, 
and will need to be taken into account explicitly. 

Finally, in this section, we determine /3c by fitting the 
crossing points of ^l/L and xsgI^^''~^ shown in Ta- 
bic VIII to Eq. (30), assuming the values in Eq. (53) and 
(54), ijj = 0.277(8), = 4.41(19). The plot is shown in 
Fig. 14, and the resuh is /3c = 0.64805(39) [2]. Combining 
the errors gives 

/3c = 0.64805(41) ^ Tc = 1.5431(10) , (55) 
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FIG. 13: (Color online) The quotients of 9^ log , 9^ log f/4 
and dfi logxsG at the crossings of xsg/L^'^~^ for a = 0.790. 
The lines represent the best straight-line fit as function of 
, using all the data, in which uj, as well as the intercept 
Q = 2^^" , is a fit parameter. 
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FIG. 14: (Color online) Values of the crossing points 
for (,l/L and Xsg / L^"^^ , for a = 0.790, as a function of 
l/L'^^^^" where the values of uj and v are fixed at the val- 
ues given in Eqs. (53) and (54). The intercept is the critical 
couphng /3c. 



with x^/dof = 4.47/10. Note that the contribution to 
the error from the uncertainty in ui is very small. 



FIG. 15: (Color online) Correlation length in units of the sys- 
tem size (top) and scale-invariant combination of the SG sus- 
ceptibility and the lattice dimension XSG / L^'^^^ (bottom), 
as a function of the inverse temperature P, for the LR-model 
with a — 0.896. For both quantities, the curves for the difi'er- 
ent L should cross at temperatures that approach the critical 
point when L grows, see Eq. (30). 



C. One-dimensional long range model with 
a = 0.896 

According to Eq. (8) and the value of 77 for the 3-d 
model given in Ref. 18, r?sR(3) = -0.375(10), a = 0.896 
is a proxy for S-d, at least according to the comparison 
of the exponents 77 (or equivalently of the magnetic expo- 
nents uh)- We now attempt to see if the correspondence 
also works for the exponents w and v. 

As we show in Fig. 15, ^l/L displays a rather marginal 
behavior for this value of a. We are not able to resolve 
the crossing temperatures for this dimensionless quantity. 
On the other hand, crossing points of Xsg/L^'^~'^ are 
easily identified. Our interpretation of these findings is 
that, for this value of cr, we are fairly close to the critical 
value <7i, such that for a > ai there is no longer a SG 
phase, see Sec. I. It is expected that^^ ai = 1 since this 
corresponds to d—2-\-r] = with d = 1 and 77 = ijLRio') = 
3 — 2(7. Hence a transition is expected for cr = 0.896. It is 
easier to find crossing points from xsg/L'^'^^^ , because, 
in the SG phase, it scales as L° with an exponent a larger 
than the corresponding one for ^l/L, so we feel that our 
results for cr = 0.896 are consistent with the expected 
transition. 
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FIG. 16: (Color online) Quotients of the dimensionless quan- 
tities ^l/L,U4 and U22 at the crossings of xsg/L^"^'^ for 
a = 0.896. Compared to the error bars there is very little 
size dependence so the data is inadequate to determine the 
correction to scaling exponent uj. 

Unfortunately, plots of dimensionless quantities do not 
allow us to determine ui because there is very little size 
dependence in the quotients. This is illustrated in Fig. 16 
which shows quotients of ^l/L, C/4 and U22 at crossings 
oixsG/L^''-^. 

To determine v we first consider the quotients of the 
logarithmic derivatives with respect to /? of the dimen- 
sionless quantities and U4 at the xsg/L^'^~^ crossing. 
A fit to Q + BiL^", does not allow us to find to, so we 
fix the value w = 0.33(3), obtained from Eq. (5) and 
the result of Hasenbusch ct al.^^ that losrI'^) ~ 1-0(1), 
obtaining 

g = 2^/" = 1.0890(202) [2] , xVdof = 1.14/5 , (56) 

which determines u to be in the range 5.7 < v < 10.4. 
Notice the smallness of the error bars coming from the uj 
error, or conversely the difficulty of determining uj from 
these quantities. 

We also tried a more complex fit including the quo- 
tients of logarithmic derivatives of the scale invariant 
quantity xsg/ . As discussed in Sec. Ill, this deriva- 
tive (but only this one) suffers from additional scal- 
ing corrections of order i"^/". Note that, according 
to Eqs. (9) and (10) and the SR valucs^^ ojsniS) ^ 
1.0(1), i'5_r(3) — 2.45(15), we expect ojlr ~ 0-33 and 
l/i'Li?, ~ 0.14, so the corrections of order L^^/'^ are dom- 
inant. We therefore fit the data for the quotients of the 
logarithmic derivative of xsg/L^'^^^ to Eq. (42), while 
for the quotients of the logarithmic derivatives of C/4 and 
^l/L wc use Eq. (36) with yo — l/i^, which corresponds 
to B2 = 0. 

To obtain a reliable fit, wc need fix the value of uj and. 



FIG. 17: (Color online) The quotients of 9/3 log ^l, 9/3 log (74 
and dplogxsG at the crossings of xsg/ L"^"'^ for a = 0.896. 
For dp log and dp log Ui the lines are fits to functions of 
type Q + BiL'", where Q = 2^/", see Eq. (36). For xsG 
we need to consider also an L"^^" term, see Eq. (42). The cj 
value is fixed to the value 0.33 expected from the 3-d data of 
Hasenbusch et al.^* who find tJsfl(3) = 1.0(1), and Eq. (10), 
while the value oi v is a. fit parameter. 

as above, we take this to be a; = 0.33(3), obtaining 

Q = 2^/" = 1.087(199)[3], xVdof = 1.54/7 , (57) 

which determines v to be in the range 6.8 < v < 10.6, so 
our estimate for v is 

j/Li?(0.896) = 8.7(1.9). (58) 

Again, the effect of the ui uncertainty is very small. We 
have tried to bound the uj value from this fit, but the 
result is almost useless [uj G (0,0.97)]. 

Finally we discuss the value of /3c. We do not see any 
evolution of /3c with L. However wc perform several fits 
to estimate the extrapolation errors. First we try a fit 
of the Xsg / L'^'^^^ crossings taking uj and v from the i-d 
derived values: uj = 0.33(3), v = 7.35(45) so uj + l/v = 
0.47(4). The result is /3c = 1.004(15)[1], with xVdof = 
0.24/2. If we use uj = 0.33(3) but the v value obtained 
above v = 8.7(1.9), i.e. uj + \/v ^ 0.44(5) we get /3c = 
1.003(16)[2],xVdof= 0.23/2. 

These last two results are statistically correlated, and 
we take the latter as our final estimate: 

/3c = 1.003(18) ^ Tc = 0.997(18) . (59) 

VI. CONCLUSIONS 

The purpose of this paper is to see if there is a value of 
a for the long-range spin glass model which corresponds 
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TABLE IX: Summary of results for critical exponents of the 
short-range models in 3-d and 4-d, the expected (proxy) re- 
sults for the long-range models based on the short-range re- 
sults and the connection in Eq. (5), and the actual results for 
the long-range models. It was not possible to estimate u) for 
the long-range model with a = 0.896. If we assume that it 
is given by the matching formula, uisB./d, then we obtain the 
result for !/lr (0.896) shown in the table. The 3-d results are 
from Ref 18, and all other results are from the present work. 





d = 4,(7 = 0.790 


d = 3,(7 = 0.896 


'^SR(d) 


1.04(10) 


1.0(1) 


LOSK{d)/d 


0.26(4) 


0.33(3) 


t»JLR(cr) 


0.277(8) 




!^SR(d) 


1.068(5) 


2.45(15) 


di^SR(d) 


4.272(20) 


7.35(45) 


^'lr(o-) 


4.41(19) 


8.7(1.9) 



nection, since, for the corresponding long-range model, 
a = 0.896, we can not determine w. If we assume 
that the value of ojlr (0.896) is that given by the match- 
ing formula, Eq. (10), with the value of uj from the 3- 
d simulations,^® namely ajLR,(0.896) = 0.33(3), then we 
find i^LR = 8.7 ± 1.9 which is consistent with 3i'sr(3) = 
7.35 ±0.45. 

While it seems unlikely to us that all the critical expo- 
nents of the LR and SR models match exactly according 
to Eq. (5), our results indicate that these equations are 
satisfied to a good approximation, and hence the crit- 
ical behavior of the SR and corresponding LR models 
are very similar. Whether this similarity extends to the 
more subtle question of the nature of the spin glass phase 
below Tc remains to be seen. 



precisely to a short-range four-dimensional spin glass, 
and (with a different value of a) to a three-dimensional 
spin glass, in the sense that all the LR and SR exponents, 
in particular, r/, v and w, match in the sense of Eqs. (5)~ 
(10). Since ?7lr is given exactly by the simple expression 
in Eq. (7), we have chosen two values of cr, 0.790 and 
0.896, as proxies for A-d and 3-c? respectively, since the 
values of 77 match according to Eq. (8). The question, 
then, is whether the other exponents, u and t^, match 
according to Eqs. (10) and (9). 

Our results for oj and v are summarized in Table IX. 
For the case of 4-d, the correspondence works well, the 
values for the exponents being consistent with Eqs. (9) 
and (10) within reasonably modest error bars. How- 
ever, for 2)-d, we are not able to establish a sharp con- 
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